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Abstract 

We have made a calculation of the flux of secondary antiprotons expected 
for the leaKyobox model and for the closed galaxy model of Peters and Mestergaard 
(1977). The p/p ratio observed at several GeV Is a factor of 4 higher than the 
prediction for the leaky-box model but Is consistent with that predicted for the 
closed galaxy model. New low energy data Is not consistent with either model. 

The possibility of a primary antiproton component Is discussed. 
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1. Introduction 

Irrespectivo of whither or not there are primary antiprotons In the 
cosmic rays* a secondary component resulting from the nuclear Interaction 
of cosmic rays In the Interstellar medium will exist. The magnitude of 
this secondary component depends critically on how cosmic rays are stored 
!n and propagate through the galaxy. 

The first reliable predictions of the anti proton to proton ratio (p/p) 
were made by Galsser and Maurer (1973) shortly after results on p production 
rates In pp collisions became available from the CERN Intersecting storage 
rings (for referenceSt see Galsser and Maurer, 1973). Their prediction was 
for the 'homogeneous* or *1eaky>box' model of cosmic ray propagation In which the 
rates of production and probabilities of escape of cosmic rays from the galaxy 
are uniform throughout the containment volume. In this model, the p/p ratio 
depends only on the shape of the proton spectrum and on the mean escape 
length, Xg. The prediction of Galsser and Maurer (1973) for an proton 
spectrum (E Is total energy) and a mean escape length of 5 g/cm (the Inter- 
stellar medium was taken to be composed entirely of hydrogen) Is shown In Figure 1. 
Later predictions by Badhwar et a1. (1975) appear to be In error as pointed out by 

Szabelski, Mdowczyk and Wolfendale (1980) and Tan and Ng (1981). 

Steigman (1977) emphasized the astrophysical and cosmological significance 
of cosmic ray antiprotons and made estimates of the p/p ratio expected for 
the closed galaxy models of Rasmussen and Peters (1975) and Peters and 
Uestergaard (1977) using the (Incorrect?) Badhwar et a1 . (1975) predictions. 

The predicted p/p ratio varied by up to a factor of 10 between the leaky-box 
and closed galaxy models. 

The first cosmic ray antiproton observations were made by Golden et a1 . 

(1979) and by Bogomolov et a1. (1979) and were found to give a p/p ratio 
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significantly higher than that predicted for the 1eaky«box model. The recent 
low energy result of Buffington, Schindler and Pennypacker (1981) confirms 
this trend. These data are shown In Figure 1. Szabelski et a1. (1980) have 
used the observed high energy p/p ratio In conjunction with their calculations 
to estimate the mean escape length of anti protons for the 1eaky*box model. 

They found this to be a factor of 5 to 10 higher than that Inferred from the 
ratio of light to medium cosmic ray nuclei (eg. boron/carbon). From this 
Inconsistency, they suggested that the distribution of pathlengths may not 
be exponential (1eaky*box model), but might have a longer 'tali' . 

Galsser, Owens and Steigman (1981) suggested that a re-examination of the 
closed galaxy model of Peters and Westergaard (1977) was merited by the observa- 
tions. A similar conclusion was also reached by Stephens (1981a) who found that 
the observed p/p ratio was about halfWay between that predicted for the leaky-box 
model and that predicted for the Rasmussen and Peters (1975) closed galaxy 
model (Stephens, 1918b). Stephens (1981a) also considered the leaky box model 
and the nested leaky box model of Cowsik and Wilson (1973). The ratio predicted 
for the latter model was lower than for the leaky box model and gave a worse 
fit to the data. In the present paper, we consider cosmic ray antiprotons In 
the context of the closed galaxy model of Peters and Westergaard (1977). Cosmic 
ray positrons In this model are discussed In a separate paper (Protheroe, 1981). 

2. The Model 

In the version of closed galaxy model proposed by Peters and Westergaard 
(1977), hereafter PW, the sources of cosmic rays are located In the spiral arms 
of the galaxy (hereafter called region S) from which they can slowly leak out Into 
an outer containment volume (hereafter called region H) which comprises part of 
the disk, and surrouding halo. The outer boundary of region H constitutes a 
closed box from which cosmic rays can not escape. Depletion of cosmic rays 


3 


in rvglon H is dut toltiy to nucloar interactions and energy losses. Region 
H thus contains an 'old component' and Region S contain a 'young component' of 
cosmic rays. The old component also permeates region S. A schematic picture 
of the closed galaxy model is shown inset in Figure 2. 

The density of the young component of protons, (cm*^ GeV*^), is 
simply: 

''p-SV*, 

where is the rate of product .on of cosmic ray protons integrated over 
the galaxy (GaV'^s"^), t^ is the residence time (s) of protons in region 

S and is equal to the mean escape time for protons leaking from region 
S into region H. is the total voline (cm^) occupied by region S. 

The density of the old component of protons, npt is given by: 

"p»«^VV„ (2) 

where fSp is the rate of production of cosmic ray nucleons (bound and un- 
bound) integrated over the galaxy (the assumption is that most nuclei leaking 
out of region S will fragment into nucleons in region H). Vie will use f « 

(1 4 S^/Sp) « I.IS. t^^ is the mean attenuation time of protons due to nuclear 

interactions in region H, and is the total volume occupied by region H. 
Equation (2) is only valid if t^^ is significantly less than the age of the 
galaxy and If Sp is independent of time as will be discussed later. 

Ue can define an escape length x^ (g/cm } and an attenuation length 
x^t ( 9 /cm ) as the amount of matter traversed by protons in region S in time 
tj and by protons in region H in time t^^ respectively. Hence, we may obtain 
the ratio of the old to young components: 

"A 

where K is the ratio of the total mass of interstellar gas in the galaxy 
to that in region S. 
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The proton energy. E, remaining after suffering a nuclear Interaction 
Is on average about half Its Initial energy E'. The distribution In energy 
may be approximated by (eg. (kisser. 1974): 

N(E)dE % dE/E' . (4) 

For a power law cosmic ray energy spectrum with differential exponent 
-Y* Md then obtain: 

*.t “ *1 

Witrt Xf Is the nean Interectlon 1en9th, tn Interstellar 

medium comprising 90S hydrogen and lOS helium nuclei by number, the mean 
mass, <m>. Is about 2.2x10*^%, and mean Inelastic cross section, 
about 3.5x10*^cm^. Thus fdr y 2.5 (as seen from equation 5 above, the 
exact value Is not critical), we obtain x^^ n% 105 g/cm of Interstellar material. 

PW found that the observed ratio of light to medium nuclei would be con* 
si stent with their model for values of K ranging from 50 to 500 and for an 
energy dependent escape length given approximately by: 

Xj •* 15 g/cm^ (6) 

for rigidities, p, greater than 4 GV/c. Since the composition of cosmic 
rays heavier than protons Is dominated by the young component (most old 
nuclei will have fragmented Into nucleons in region H), we note that Is 
essentially Identical to the mean escape length, x^, derived using the leaky* 
box model. We therefore use the result of a more recent analysis by 
Protheroe, Ormas and Comstock (1981) using the 1eaky*box model: 

7 g/OB^^ pj 4 Gv/c . . 

^s •* 7(J)“°*^ g/cm^, p > 4 GV/c. 

We have used equations (3), (5) and (7) to decompose the observed proton 
spectrum Into Its young and old components. For the total proton spectrum, we 
have used the range of demodulated spectra from the work of Morflll, VDik and 
Lee (1976). This range of spectra is plotted as the shaded area in figure 2. 
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Ue also shoM the old component obtained as described above for values 
of K ranging from SO to 500 and the young component for K > 100. 


3. Anti proton Production 

Having obtained the cosmic ray proton spectra In the tMO regions of 

the closed galaxy model, we proceed to calculate the rates of production 

of anti protons In these two regions. First we must consider the production 

of antiprotons In high energy pp Interactions. 

It Is convenient to describe the Inclusive reaction, p p p 

'anything'. In terms of the CM energy stquared. s. and the Fbynman scaling 

variable, x • 2p»*//s“ (* refers to the CM system). Then, following Galsser 

and Maurer (1973) we define a function F - In terms of the Invariant cross 

PP 

section for the process, E ^7 

f- (x,s) . ^iST“ 2.P, E 0 dp^. (8) 

rr «pp 0 V u 

where Is the maximum value of the transverse momentum, p^, which is 
kinematically allowed for given values of x and s. 

In the most recent results from Femrilab (Johnson et a1., 1978), the 
Invariant cross section at fixed values of p^ and the radial scaling variable, 
X|^, was found to be Independent of s over the range of incident proton 
energies available at Femrilab. The radial scaling variable Is defined by: 


where Is the maximum energy In the CM system that Is available to an 
Individual antiproton and Is given by (Taylor et al., 1976): 


(9) 


Ei., •• (s - 7.02)/2/? GeV. (10) 

ihBa 

Xp Is thus a function of x, s and p^, except that Xj^-* |x| as s 

The data of Johnson et al. (1978) are well fitted by a function of the 
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fbrn: 

^ (l+P^Vb)** 

with; A • 1.9 +0.4 rtb/teV^; n ■ 8.1 il.4; and b • 1.2 ^0.3 (GeV/c)^. In 
Figure 3 we plot the product, Opp^"*^(s) Fp^x.s), obtained from these 
fits as a function of x for various values of s. Ue point out that the 
dependence on s at x • 0 extrapolated from the Ferarilab data Is consistent 
with that observed over the ISR energy range by Guettler et a1. (1976) as 


( 11 ) 


well as with the low energy data as fitted by Galsser and Maurer (1973), 

These ISR data may therefore not Indicate a breakdown of scaling In the 
central region (|x| « 1), but Instead that the Feynman scaling limit may 
not yet have been reached (see also Ellsworth, 1979). 

The rate of production of anti protons of energy E, per Interstellar 
nucleon, P(E), In cosmic ray nuclear Interactions In the interstellar medlun 
Is given by: 


P(E) • 2 r" M(r)Np(n aJj«^(s')Fp^x,s') ^deOerU-^ (12) 

^ nrin 

Here, Np is the density of cosmic ray protons (cm’^ GeV"^ and is either Hp 
(region H), or (Yp'^^p) (region S); is the minimum proton energy required 
to produce an anti proton of energy E; s' is the square of the CM system 
energy, 2n^(P * n^); and M is a correction factor to take into account p*He 
interactions (regions S and H), and at>H interactions (region S only). This 
factor has been discussed by Giler et a1. (1977), and we use here. values 
applicable to charged pion production from our earlier work (Protheroe, 1981b). 
The factor of 2 in equation (12) is to include the production of antineutrons 
which decay into antiprotons. We approximate x in equation (12) by (Gaisser 
and Maurer, 1973): 


X m [(£^-rap^-<P^^>)** - E(l-4 mp^/s)'*]/mp 


( 13 ) 


where the mean square transverse momentum was obtaine.'i from the Johnson 
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et a1. (1978) fits: 

1-0.59 (6tV/c)^. (14) 

The rates of production of anti protons obtained as described above for the 
tMO regions of the closed galaxy model are plotted In Figure 4. Fbr region 
H, we show the production rates expected for values of K ranging from SO 
to 500. 

For comparison with previous work, we have calculated tha p/p ratio 

2 

for the leaky-box model with an escape length of 5 g/cm of hydrogen and 

-2 6 

for a proton spectrum of the form E : 

p/p L p > W”'** ^ 

^ ^nrin 


This Is plotted In Figure 1 and found tc be In good agreement with recent 
results (Szabelski et a1., 1980; Tan and Ng, 1981). 

4. Eoulllbrium Antlproton Spectrum 

To obtain the equillbriun spectrum of antiprotons one must take Into 
account. In addition to their rates of production and escape from the 
containment volume, the effects of Ionization losses and energy losses due 
to nuclear Interactions (equation 4). The equilibrium density of antiprotons 
of energy E, Np<E), Is then obtained from: 

d^ 




(16) 


where: x^ Is the mean escape length (x^ for region S; • for region H); 

and X. are the mean lengths for nuclear Interaction and annihilation; and 
on 

2 

dE/dx Is the rate of ionization energy loss (GeV per g/cm }. While measure* 

ments of the total Inelastic cross section, , have been made up to 

PP 

100 GeV, the annihilation cross section has only been measured up to kinetic 
energies of 'v< 6 GeV (Spencer and Edwards, 1970; Locken and Derrick, 1963; 
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Anidi ttt a1., 1966; Lynch et a1., 1963; Brou, 1973; Bockmann at a1.. 1966; 

Fartcl et a1., 1968). To obtain the annihilation and non-annihilation In- 
elastic cross sections* and o-p^* separately at higher energies, m shall 
adopt tMO different approaches. The first Is to extrapolate the observed 
annihilation cross section to high energies using a 1/$* dependence (black 
absorbing sphere eodel: Kobe and Takeda, 1958). This gives »> 24/(1-»‘2iyT)%t. 
The measured as_ Is tiien consistent with o- * (o-_ ♦ os- ^) with 

pp pp ' pp pp ' 

os- ^ 10 mb. The second approach Is to assume o- ^ and to obtain 

pp PP PP 

o-p*"by subtracting o-p ^ from o-p This leads for the energy range of 

Interest to os. ^ * (32 - 5T “**)mb and o-_ w 36T"** mb. This latter fona for 
the annihilation cross section Is also consistent with the available data and mav 


be more consistent with expectations of quark models at high energies 

(T. K. Galsser, private connunl cation). To obtain the annihilation and Inter- 


action lengths, x « <m>/os , we assume o- * 1.07 o- . The term 

P - ISM P - ISM PP 


Involving the Integral In equation (16) Is to account for Inelastic Interactions 


In which the p does not annihilate. 


It should be pointed out that the steady state solutions used here 
(equations 2 and 16) are only valid If the attentuatlon time of protons 
In region H, t^^, and the Interaction-annihilation time of anti protons, 

(1/x^ ♦ '*** 9P^«k>'» 

Also, the conditions In the galaxy (matter density, ccssilc ray sources, etc.) 
should have changed little during this time. The fl^t conditions are not 
mat If the mean nunber density of Interstellar nucleons, n^^. In region H Is 
less than -v 3x10'^ cm"^. In this case, the predicted anti proton flux will 
be too high and should be taken as an upper limit. 
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Equation (16) has boon solvad for rt 9 lon S to ylald the young ooiqMnont 

« 

of tha antiproton flux uhlch Is plottad In Pigura 5. Nota that this flux 
Is Idantical to that pradictad for tha laaky-box nodal with naan ascapa 
langth givan by aquation (7). Fbr raglon H, tha antiproton production rata* 
P* In aquation (16) contains an additional tans to account for laakaga fron 

a 

raglon S Into raglon H. This additional sourca Is givan by: 




and Is plottad In Figure 4 for tha valuas of K considarad bars. Equation 

(16) has than baan solvad for raglon H to ylald tha old conponant of anti- 

proton flux and this Is plottad In Figure 5 for values of K ranging frosi SO 

to SOO. tfa show tha affict hare of using diffarent extrapolations of 
on tha predicted flux for K ■ SO. 

Since ttw young component of tha antiproton flux Is the sama fOr all valuas 
of K (It Is datemrinad only by tha proton spectrum In tha spiral arms and tha 
ascapa length), tha variation of total p flux with K Is datarmlnad by tha old 
component. Frwn figure 4 wa sea that tha production rate par nucleon Is higher 
fbr low values of K. This Is expected from our decomposition of the proton 
spectrum Into Its young and old conponents (Figure 2) described by equation (3). 
Thus In figure S wa sea also that for high values of K the old conponant of 
antiproton flux Is lowest. 

S. Conclusion 

Tha p/p ratio obtained by dividing the sum cf the young and old components 
of tha antiproton spectrum (given In Figure S) by the 'demodulated proton spec- 
true* of Figure 2 Is plottad In Figure 6 fbr different values of K, together 
with tha prediction for tha laaky-box model. The decrease with energy at high 
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tntrglts of tht p/p ritlo for th« Itaky box nodtl Is duo to our use of an energy 
dependent escape leng^. Tht prediction of S^hms (1M1) Is also shown for 
the closed galu^ andel of Rasmussen and Peters (1975). The p flux predicted 
for this model Is equivalent to ^ g]l component of the flux In tht PW model 
for K • 1 lit Is Interesting to note that our leaky box model prediction Is 
equivalent to In the PU model). Stephens' prediction for the leaky box 
model with energy dependent escape (not shown) Is also consistent with our result. 
We compare these predictions with tiie observed ratios In the figure. The high 
energy data (Golden et a1., 1979; Bogomolov et a1.. 1979) are fbund to be 
consistent with the predictions of the closed galaxy model for the samr range 
of K that Is required to reproduce the observed ratio of light to medium nuclei 
In the msmic rays. Other aspects of this model should therefore be exaarined In 
detail (eg. positrons* Y*rays* ratio of light to medium nuclei at high energies). 
Such a study has already hiien undertaken for the case of positrons (Protheroe* 
1331). Fror this *^ork, it .nas fo-Td tiat t.n orsJictec oositron flux for the 
closed relexy model with K • lOO was consistent with the 'resen: Jata provided 
the Interstellar density In region H* n^« was less than 0.3 nucleons/cm^. 
Attempts should also be made to put the model on a more soimd theoretical 
basis by seeking a feasible mechanism to 'close the galaxy' to cosmic rays, 
at least up to 100 GeV. 

The low energy data (Buffington et al., 1981) are more difficult to 
fit by secondary antiproton production. In Figure 6 we also show the effect 
of solar modulation on the p/p ratio for K ■ SO. For this, we have assumed 
that the modulation at the near solar maxlmun conditions when the observations 
were made may be described by a mor ^mrqy loss of 600 NeV (Urch and Gleeson, 
1973). The observed ratio Is still more than a factor of 10 (3.5 standard 
.iQvlatlons) above >he predicticn. 
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An altnmatlvt txplanatlon In tarm of a primary antiproton cooponont 
(Qalsstr and Ltvy* 1974; StodMr, 1981) should not bt discounttd. If tho 
Itaky box modal with an anargy dapandant mtan atcapa langth (aquation 7) 

Is a battar dasciiptlon of tha mal situation than tha closad galaxy modal* 

•A 

than a primary antiproton c omp on ant of magnituda <>> 3x10 of tha otaaruad 
proton spactrum would axist. This possibility Is discussad further by 
Stacker* Protharoa and Kazanas (1981). 
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Floura Captions ; 

Figure 1: Predictions of p/p ratio by Galsser and Maurer (1973) and 

later authors for leaky-box nodel with mean escape length 

of 5 g/cm of hydrogen Independent of energy and for a proton 

•2 6 

spectrum proportional to E . Observations of Golden et a1. 
(1979) (f). Bogomolov et a1. (1979) (0) and Buffington et a1. 
(1981) (m) are Indicated. 

Figure 2: The observed proton spectnmi as decomposed Into Its young 

and old components for different values of K. Shown Inset 
Is a schematic representation of the regions of the galaxy 
In the closed galaxy model (not to scale). 

Figure 3; The dependence of Ppp (equation 8) on x for various 

values of s. 

Figure 4: The rates of production of anti protons In regions S and H as 

a function of kinetic energy. Note that production rate for 
region S Is Identical to that for leaky box model. The source 
term due to leakage from region S Into region H (equation 17) 
Is also shown. 

Figure 5: Predicted young and old components of antiproton flux. 

Figure 6: Predicted p/p ratio for 1eaky>box model with energy dependent 

escape and the closed galaxy model of Peters and Westergaard 
for various values of K. The dashed line is the prediction 
of Stephens (1981b) for the earlier closed galaxy model of 
Rasmussen and Peters. The effect of solar modulation c>i the 
prediction for K • 50 is shown by the curve labelled K ■ 50 
(mod.). See Figure 1 for key to data and Figure 5 for key 
to annihilation cross section used. 
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